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Abstract 

In this article, we explore the interplay between searches for supersymmetric particles and Higgs 
bosons at hadron colliders (the Tevatron and the LHC) and direct dark matter searches (such as CDMS, 
ZEPLIN, XENON, EDELWEISS, CRESST, WARP and others). We focus on cohider searches for 
heavy MSSM Higgs bosons {A, H, H^) and how the prospects for these searches are impacted by direct 
dark matter limits and vice versa. We find that the prospects of these two experimental programs 
are highly interrelated. A positive detection of A, H or at the Tevatron would dramatically 
enhance the prospects for a near future direct discovery of neutralino dark matter. Similarly, a 
positive direct detection of neutralino dark matter would enhance the prospects of discovering heavy 
MSSM Higgs bosons at the Tevatron or the LHC. Combining the information obtained from both 
types of experimental searches will enable us to learn more about the nature of supersymmetry. 
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I. INTRODUCTION 



Recently, direct dark matter detection experiments, most notably CDMS, have reached 
the level of sensitivity needed to detect neutralino dark matter over a substantial range of 
supersymmetric parameter space 0, 0, 0, lEI . The reach of these experiments continues to 
expand rapidly, probing an ever increasing fraction of supersymmetric models. 



As searches for neutralino dark matter progress pi], searches for supersymmetry at Run II 
of the Tevatron are becoming increasingly sensitive as greater integrated luminosity is accumu- 
lated. Results from Tevatron searches for neutralinos and charginos 0], MSSM Higgs bosons 
[sl and squarks and gluinos [§] have recently been published. Although no evidence for su- 
persymmetry has yet been found at the Tevatron, many of these limits are the strongest to 
date and continue to advance. Furthermore, the first collisions at the Large Hadron Collider 
(LHC), with a center-of-mass energy of 14 TeV, are scheduled to take place in 2008. If low 
energy supersymmetry exists in nature, it will likely be discovered at these experiments in the 
relatively near future. 

These two very different experimental programs - dark matter searches and collider exper- 
iments - are carried out largely independently of each other, with relatively little interaction 
between their respective communities. There is, however, a great deal of physics interplay 
between these two search strategies. Direct dark matter detection prospects depend on the 
mass and couplings of the lightest neutralino, as well as those of the Higgs bosons and squarks 
exchanged in elastic scattering diagrams. Many of these properties can potentially be measured 
or constrained by collider experiments. Turning this around, an astrophysical detection of neu- 
tralino dark matter, or an upper limit on its elastic scattering cross section with nucleons, can 
be used to provide information which is valuable to the accelerator community. 

In this paper, we expand on earlier work |10|, which explored the interplay between direct 
dark matter experiments and searches for supersymmetric Higgs bosons at the Tevatron. In 
this article, we focus on Tevatron and LHC searches for heavy MSSM Higgs bosons {A, H, H^), 
and discuss their relationship to direct neutralino dark matter searches. The paper is organized 
as follows. In Sec. [TTlwe introduce the aspects of the supersymmetric spectrum relevant to our 
study. In Sec. Illll we review the contributions to the spin-independent neutralino-nucleon elastic 
scattering cross section and explore the impact on this cross section of the supersymmetric 
particle spectrum relevant to collider searches. In Sec.[IV]we present the current exclusion limits 
and future discovery reach of the Tevatron and the LHC in the search for heavy MSSM neutral 
Higgs bosons, and compare this to the reach of direct neutralino dark matter searches. In Sec. |V] 
we discuss the implications and interplay of these two very different classes of experiments. In 
Sec. I VI I we describe the caveats and limitations of our analysis. Finally, we summarize and 
state our conclusions in Sec. IVIII 

II. NEUTRALINOS, SQUARKS AND HIGGS BOSONS IN THE MSSM 

In the Minimal Supersymmetric Standard Model (MSSM), the superpartners of the four 
Standard Model neutral bosons (the bino, wino and two neutral higgsinos) mix into four physical 
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states known as neutralinos. The neutralino mass matrix in the B-W-H1-H2 basis is given by 



Mn = 



/Ml -Mz cos /3 sin 6'y[/ Mz sin (3 sin 6^, \ 

M2 Mz cos P cos 9^^ — Mz sin P cos 6^r 

—Mz cos P sin Mz cos P cos 9^/ — 

^ Mz sin P sin —Mz sin P cos9^r -/jl y 



(1) 



where Mi and M2 are the bino and wino masses, n is the higgsino mass parameter, 9w is the 
Weinberg angle, and tan P is the ratio of the vacuum expectation values of the up and down 
Higgs doublets. This matrix can be diagonalized into mass eigenstates by 

^diag ^ N^M^oN. (2) 

We are interested here in the lightest neutralino, which in the presence of R-parity conservation 
can serve as a viable dark matter candidate.^ In terms of the elements of the matrix, A^", the 
lightest neutralino is given by the following mixture of gaugino and higgsino components: 

X° = NuB + Ni2W^ + iVia^i + iVi4^2. (3) 

The mass and composition of the lightest neutralino is a function of four supersymmetric 
parameters: Mi, A'/2- /i and tan/5. This becomes further simplified if the gaugino masses are 
assumed to evolve to a single value at the GUT scale, yielding a ratio at the electroweak scale 
of Ml = I tan^ 9wM2 ^ O.5M2. 

As we will describe in the following section, neutralinos can scatter with nuclei through the s- 
channel exchange of the superpartners of the Standard Model quarks. These contribution can be 
particularly important in the case of light squark masses and large to moderate values of tan/?. 
The diagonal entries in the squark mass matrices are largely set by the soft supersymmetry 
breaking mass parameters mQ^, mu^ and mo^, where i — 1, 2, 3 denotes the generation. In many 
supersymmetry breaking scenarios, such as minimal supergravity and gauge mediated models, 
these parameters are nearly degenerate, leading to all of the first and second generation squarks 
to naturally have very similar masses. Third generation squarks (stops and sbottoms) have their 
masses split by their Yukawa couplings, however. The stop and sbottom mass matrices are given 
by: 

Ay/2 ^ f^Q3+^t+ ^UL (A - cot P) \ , . 

\mt {At - II cot P) mfj^ + mf + D^^ J 



and 



^ ( "^^3 +^b + rub {Ab - /i tan P) \ 
b \mb {Ab - tan P) m^.^ + ml + D^^ J ^ ' 

Here, the D's are D-term contributions, and are of order m|. The quantities At and Ab are soft 
supersymmetry breaking trilinear scalar couplings in one-to-one correspondence to the Yukawa 
terms in the superpotential. The mass splitting in the stop sector can be especially pronounced. 
Defining a common squark mass rriq = mg. = m[/. — niDi, then for rriq ^ rrit, the physical stop 
masses are well approximated by m^^ ^ ^ [m~ ^ mt{At — ficot P)]^^"^. 



^ R-parity is defined as i? = (— i)3S+l+2S^ where B, L and S denote baryon number, lepton number and spin, 
respectively, i? = 1 for all SM particles and -1 for their superpartners. 
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Throughout this study, we will consider two benchmark scenarios which provide conventional 



choices for the mixing in the stop sector The first of these, called the "no-mixing scenario" , 
is defined hj At = n cot f3, minimizing the mixing in the stop sector and yielding the minimal 
radiative corrections to the lightest CP-even Higgs mass, ruh. The second case, called the "m™*^^ 
scenario", is defined by At = 2m^ + /icot (3, and maximizes the value of m/j. For definiteness, 
we adopt A^ = At throughout. 

In addition to exchanging squarks, the elastic scattering of neutralinos with nuclei can be 
induced through the exchange of Higgs bosons. Rather than the single Higgs boson of the 
Standard Model, a minimum of two Higgs doublets are required in supersymmetric models 
to avoid triangle diagram anomalies and to create a gauge-invariant superpotential that gives 
masses to both up and down-type quark fermions. In the absence of explicit CP-violation, 
the two Higgs doublets of the MSSM correspond to physical states in the form of two neutral, 
CP-even Higgs bosons [h and H), one neutral CP-odd Higgs boson [A) and one charged Higgs 
boson (H^). The angle a diagonalizes the CP-even Higgs squared mass matrix. 

It is traditional to take as a free parameter which, together with the value of tan/3, 
determines at tree level the masses of all the other Higgs bosons: 

m^± = m\ + m^ (6) 

^Ih = ^ (^^A + \l{m\ + miy - Amlm\ cos'^ 2/?^ (7) 

In the case that rriA ^ rnz, it follows that rriH ~ itt-h^ ~ and cos a ~ 1. These simple rela- 
tionships are significantly modified at the loop level, however, and become strongly dependent 
on the supersymmetric parameters in the stop and sbottom sectors. In our study, we have used 
the Feynhiggs package |12| to calculate the Higgs masses.^ For a more complete discussion of 
Higgs and supersymmetry phenomenology, see Refs. 14 1, 15 1 and [l6|. 



III. DIRECT DETECTION OF NEUTRALINO DARK MATTER 

Experiments such as CDMS [1], ZEPLIN 0, EDELWEISS Q, CRESST ^] and WARP Q 
attempt to detect dark matter particles through their elastic scattering with nuclei. This class 
of techniques are collectively referred to as direct detection, in contrast to indirect detection 
efforts which attempt to observe the annihilation products of dark matter particles. 

Neutralinos can scatter with nuclei through both scalar (spin-independent) and axial- vector 
(spin-dependent) interactions. The experimental sensitivity to scalar couplings benefits from 
coherent scattering, which leads to cross sections and rates proportional to the square of the 
atomic mass of the target nuclei. The cross sections for axial-vector elastic scattering are 
proportional to J(J -|- 1), however, and thus do not benfit from large target nuclei. As a 
result, the current experimental sensitivity to axial- vector couplings is far below that of scalar 
interactions, and well below the range predicted for neutralinos. For this reason, we consider 
only the scalar case here. 



^ Similar results can be obtained using the CPsuperH package 13 1. 
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The scalar neutrahno-nuclei elastic scattering cross section is given by: 

■[Zfp + {A-Z)U 



a 



X 



/ nj ) 



7r(m^o + ttit) 

where is the target nuclei's mass, and Z and A are the atomic number and atomic mass of 
the nucleus, fp and /„ are the neutralino couplings to protons and neutrons, given by: 



rrir 



m, 



,71 r(p,n) \ ^ ''''p,n 

^ q=c,b,t 



q=u,d,s ^ 

where the neutralino-quark couplings and, being conservative, / 



(9) 

0.020 ± 0.004, 



/jf^ 0.026 ± 0.005, 4f 



0.118 ± 0.062, f. 



0.014 ± 0.003, /, 



(n) 



0.036 ± 0.008 and 



(n) 



0.118 ±0.062 n 



The first term in Eq. [9] corresponds to interactions with the quarks in the target nuclei, 
either through t-channel CP-even Higgs exchange, or s-channel squark exchange: 



H,h 






The second term corresponds to interactions with the gluons in the target through a 
quark/squark loop diagram, f^^ is given by 1 — f!^J — f^^ — f!j?J ^ 0.84, and analogously, 

/j!^ ^ 0.83. To account for finite momentum transfer, the calculation should also include the 
appropriate form factor. 

The neutralino-quark coupling, in which all of the SUSY model-dependent information is 
contained, is given by fl8| : 

^ -Re [(X,) (Yi)*] - -—^ -Re [{W.^ 



2(m2. -m2)' 



2(m2. - 



AmwB 



Re {5i[g2Ni2-9iNii])DC 



1 



1 



2 ~^ 2 



1)2 (J2 

+ Re {62[g2Nu - 9iNu]) ( — + — 



(10) 



where 



Vi 



, g2mqN*^_^ ^ 
* 92mgNi^5 



2mwB 



V22 



2m^B +^Hf^^^n,+..T.iV.), 



(11) 



5 



where throughout i = 1 for up-type quarks and i = 2 for down type quarks. mii,m2i denote 
the squark mass eigenvalues and r] is the matrix which diagonahzes the squark mass matrices, 
diag{m\,m\) = rjM'^ri'^. i/i, T^i and Cj denote hypercharge, isospin and electric charge of the 
quarks. For scattering off of up-type quarks: 

5i = A''i3, 62 = Ni4, 5 = sin/3, C = sina, D = cosa, (12) 

whereas for down-type quarks: 

Si = Ni4, 62 = —Ni3, 5 = cos/5, C = cosa, D = —sma. (13) 

The first two terms of Eq. [10] correspond to interactions through the exchange of a squark, 
while the final term is generated through Higgs exchange. 

We will now describe the behavior of this cross section in a number of interesting limiting 
cases. Firstly, in the case of heavy squarks, small wino component (A'^12) and cosa ~ 1 (which 
implies moderate to large tan/5 and rriA ~ mn ~ ^h±)i elastic scattering is dominanted by H 
exchange with strange and bottom quarks, leading to a neutralino-nucleon cross section of: 

^ 4vrm^cos^/?m|, 1,^^^ + 27^^^ ) ' ^"^^ ^"'^e, cos a - 1). (14) 

If instead, we consider the case of heavy squarks, small A^i2 and heavy H (cosa ~ 1), we 
find that scattering with up-type quarks through light Higgs exchange dominates: 

(^xN 47rm^ V 27 J ' " ~ y^^) 

Thirdly, for the case of scattering dominanted by squark exchange, with large to moderate 
tan (3 and approximately diagonal squark mass matrices, we find the limit of: 

We can see these behaviors exhibited in Fig. [1] where we plot the neutralino-nucleon spin- 
independent elastic scattering cross section as a function of the CP-odd Higgs mass, (recall, 
that for niA ^ ~ ^h)- In this figure we have adopted the no-mixing scenario and 

consider different values of tan/5 and ^.^ For relatively small values of m^, the cross section 
scales with l/m\, as found in the limit of Eq. [TH As itia becomes larger, the cross section 
becomes dominanted by h exchange and flattens with respect to m^, as found in the limit of 
Eq.[l5l 

For large tan /5, large itia and light squark masses (the limit of Eq. [T6|) , the effect of the squark 
exchange contribution on the elastic scattering cross section becomes significant. In Fig. [2] we 
plot the neutralino-nucleon spin-independent elastic scattering cross section for various values 
of tan (5 and /x as a function of the squark masses. From this we see that in the case that the 
H exchange contribution is suppressed (due to large ^ ""^h), squark exchange diagrams 
dominate the elastic scattering cross section if the squarks are light and tan (3 is large. 



^ Note that in the case of the no- mixing scenario, the tan/3 = 3 contours appearing in Figs. [T] and [2] predict a 
value of nih which is excluded by LEP. 
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FIG. 1: The neutralino-proton spin- independent elastic scattering cross section as a function of 
for selected values of n (200, 400, 800 and 1200 GeV) and tan/? (3, 10, 30 and 60). For each case, 
Ml = 100 GeV, M2 = 200 GeV, At = Ab = fi cot /3 (the no-mixing scenario) and = 1 TeV have 
been used. The results found for At = Ah = 2mq + //cot (3 (the m™*™ scenario) are very similar. The 
solid lines do not include any contributions from squark exchange, while the dotted lines include both 
Higgs and squark exchange. Shown as a horizontal dashed line is the current upper limit from the 
CDMS experiment [l|. 
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FIG. 2: The neutralino-proton spin-independent elastic scattering cross section as a function of rriq 
for selected values of fi (200, 400, 800 and 1200 GeV) and tan/? (3, 10, 30 and 60). For each case, 
Ml = 100 GeV, M2 = 200 GeV, At = Ab = fi cot P (the no-mixing scenario) and ttia = 1 TeV have 
been used. The results found for At = A^ = 2mq + /icot fi (the m^^^ scenario) are very similar. The 
solid lines do not include any contributions from Higgs exchange, while the dotted lines include both 
Higgs and squark exchange. Shown as a horizontal dashed line is the current upper limit from the 
CDMS experiment [l|. 
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IV. HEAVY MSSM HIGGS BOSON SEARCHES AT THE TEVATRON AND LHC 



Searches for heavy MSSM Higgs bosons are being carried out at Run II of the Tevatron. In 
particular, efforts are underway to observe the processes: 

pp A/H + X ^ T+T- + X, 

pp ^ A/ H + bb ^ bb + bb, (3b's tagged) 

pp ^tt^H'^ + bb T^Ur + W"^ bb, 

pp A/H + b^ T^T- + b. (17) 

At the LHC, such particles may be observed via 

pp A/H + X ^ T+T- + X, (18) 
pp ^ H^ + tX ^ T^Ur + tX. (19) 

In each of these cases, the prospects for discovery are much greater in the case of large tan (3 
and relatively light m^. In the large tan/5 regime, the leading contributions to the production 
of an y4 or if rely on the b-quark Yukawa coupling, and thus scale with tan^ j3. For example, 
A/H production via gluon fusion is dominanted by diagrams with a b-quark loop. The cross 
section for charged Higgs production at the LHC also scales with tan^ f3. 

The current limits on heavy MSSM Higgs bosons from Run II of the Tevatron have been 
published for each of the channels shown in Eq. [T71 CDF has published limits for the di-tau and 
charged Higgs channels using their first 310 and 193 pb~^, respectively [l^. DO has published 
(or presented) limits on the inclusive di-tau, bb+bb and 6-tau-tau channels using their first 348, 
260 and 344 pb~^, respectively [23]. The Tevatron currently has approximately 1 fb~^ of data, 
and is expected to accumulate a total of 4 fb~^ or more by the end of its operation. 

The results of the heavy MSSM Higgs searches are generally represented in the tan/?-m^ 
plane, for a given choice of /j, and other SUSY parameters (which impact the limit through 
radiative corrections [2l|). The limits found for the di-tau channel, at both the Tevatron and the 
LHC, are the most interesting and are quite robust to variations in fi and other supersymmetric 
parameters. In the special case of negative and large n, the limits from the b-quark channel 



at the Tevatron can be somewhat more constraining than the di-tau channel [2l[, although 
this is disfavored by the combination of measurements of the muon's magnetic moment 
and the B —>■ Xgj branching fraction [3, 2^. The limits on charged Higgs bosons from the 
Tevatron, and as projected from the LHC, yield weaker constraints in the tan (3-mA plane at 
this time. For these reasons, we here focus on the di-tau channel for heavy Higgs searches at 
both colliders. 

Since both heavy MSSM Higgs searches and neutralino direct detection depend strongly on 
iriA and tan/3, we can compare the limits and projected reach of these experiments to each 
other in this plane [13]. In Figs. [3]and|ll we show the current and projected limits from direct 
dark matter experiments (CDMS and Super-CDMS) along side the current and projected limits 
for heavy Higgs searches at the Tevatron and LHC in the A/H + X t^t~ + X channel, for 
the m™^^ and no- mixing benchmark scenarios (see Sec. [ITl)- From these plots, it is clear that 
those regions of supersymmetric parameter space most accessible to heavy Higgs searches at 
the Tevatron and LHC are also most likely to be probed by direct detection experiments, and 
vice versa. The prospects for heavy MSSM Higgs boson detection at colliders are, therefore, 
quite correlated to the prospects of direct detection experiments. 
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As can be seen in Figs. |3] and HJ the main difference in the prospects for these two types 
of searches is that direct detection rates depend critically on the values of /i, Mi and M2 
through their impact on the composition and mass of the lightest neutralino, while the reach of 
the inclusive r+r^ channel at the Tevatron/LHC is largely independent of the precise values of 
these parameters. Smaller values of /x relative to Mi yield larger elastic scattering cross sections 
and, therefore, lead to much greater sensitivity in the m^-tan/? plane. The dependance on these 
results of the other supersymmetric parameters are generally mild. For example, the results 
shown in Fig. [3] (the m™'^^ scenario) are very similar to those found in Fig. H] (the no-mixing 
scenario) . 



10 




FIG. 3: The range of parameters in the tan (3~mA plane that can be excluded by present and future 
direct detection experiments (for various values of /x), compared to the exclusion reach for heavy 
Higgs searches at the Tevatron and LHC. Regions to the left of the contours are (or will be in the 
future) excluded by the given experiment. The upper left frame compares the current bound from 
CDMS to the 1 fb"^ exclusion region of the Tevatron. The upper right frame compares the 2007 
projected limit from CDMS to the projected Tevatron limit after 4 fb^^ integrated luminosity. The 
lower frames compare the projected limits of Super-CDMS (phase A and phase C in the left and right 
frames, respectively) to that of the LHC with 30 fb-^ In each frame, Mi = 100 GeV, M2 = 200 GeV, 
At = 2mq + n cot P (the m^^^ scenario) and niq = 1 TeV have been used. The direct detection limits 
arc shown as solid lines which do not include contributions from squark exchange and as dotted lines 
which include both Higgs and squark exchange. The values of jj, used for calculating the Tevatron/LHC 
reach are ^=200 GeV (dotted), /i=400 GeV (dot-dash), ;u=800 GeV (dot-dot-dash) and /x=1200 GeV 
(long dashed). 
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FIG. 4: The same as in Fig. [3l but with At = ficot (3 (the no- mixing scenario). 
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V. THE IMPLICATIONS OF DIRECT DETECTION FOR COLLIDERS AND VICE 
VERSA 



As collider and direct detection experiments continue to operate, their results will have 
considerable implications for each other. Consider, for example, the case in which a positive 
detection of neutralino dark matter is made in the near future by CDMS. For models with a 
sufficiently large cross section to be detected by CDMS, the neutralino-nucleon scattering is 
typically dominanted by heavy Higgs exchange, and can be reasonably well approximated by 
Eq. [m A large cross section therefore would be expected to favor a small value of rriA and, to 
a lesser extent, a large value of tan/3. 

Turning this scenario around, we can imagine a case in which a positive detection of the 
process pp ^ A/H + X ^ t~^t~ + X is made in the future at the Tevatron. Such a detection 
would imply both a large value of tan f3 and a relatively small value of m^- From Figs. [3] and HJ 
it is clear that if heavy Higgses are observed at the Tevatron, then CDMS is likely to observe 
neutralinos in the near future unless the higgsino fraction of the LSP is very small (|/i| is large). 

To study the interplay between collider and direct detection searches more systematically, 
we have performed a scan over a range of supersymmetric parameters. In particular, we have 
varied M2, fi, rriA, At, Af, and the sfermion masses up to 4 TeV (and with either sign of /i) and 
tan /? between 1 and 60. We have set the parameters Mi and M3 according to the conditions for 
gaugino mass unification. In the upper left frame of Fig. [5l we have plotted, in the tan-/3 — rriA 
plane, the models found in our scan which are not currently excluded by CDMS, but are within 
the reach of CDMS projected for 2007. Although models were found over most of the range of 
tan P and shown, the majority are concentrated at low uia and high tan P, as expected.^ 

In the upper right frame of Fig. Owe plot, in the m^- a^^ plane, those models from our scan 
which are within the 3a discovery reach of the Tevatron A/H + X t^t~ + X search (with 4 
fb~^). As expected, we find that large elastic scattering cross sections are found for this subset 
of supersymmetric models. The overwhelming majority of models we find within the reach of 
the Tevatron are either currently excluded by CDMS or will be detected by CDMS in the near 
future, or by the first phase of super-CDMS. 

In the upper left and upper right frames of Fig. El we have required that each point shown 
does not violate any direct collider constraints {ie. Higgs, chargino, slepton and squark mass 
limits) and does not predict a thermal abundance of neutralinos in excess of the dark matter 
density as measured by WMAP {VL^h'^ < 0.131) [11]. If we also require that the thermal 
density of neutralinos constitute the entire dark matter abundance (as opposed to non-thermal 
contributions being substantial), then most of the models shown in the upper frames of Fig. [5] 
must be discarded. The lower left and lower right frames of Fig. [5] are similar to the upper 
frames, but only include those supersymmetric models that produce a thermal relic abundance 
within the 3a range measured by WMAP (0.131 > VL^K^ > 0.089). 



Our scan was carried out logarithmically over each parameter, and therefore we have shown the results of 
this scan over log-scale axes in Fig. [5l in constrast to our earlier figures in the m^-tan /3 plane which were 
presented with a linear scale. 
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FIG. 5: Upper left: The values of and tan /3 for randomly selected supersymmetric models which 
are not currently excluded by CDMS, but are within the reach of the 2007 CDMS projection. The reach 
of the Tevatron and LHC heavy Higgs searches are also shown. The near future detection of neutralinos 
by CDMS enhances the prospects for heavy Higgs discovery at collliders, but does not guarantee any 
such result. Upper right: The values of and a^N for randomly selected supersymmetric models 
within the 4 fb~^ 3a discovery reach of the Tevatron heavy Higgs search. The 2007 projected reach of 
CDMS is also shown. In each frame, the points shown evade current collider constraints and provide 
a thermal abundance of neutralinos that does not exceed the measured dark matter density. Lower 
left and right frames: The same as the upper frames, but only showing those models in which the 
predicted thermal abundance of neutralinos matches the dark matter density as measured by WMAP 
(0.131 > n^h"^ > 0.089) [2^. (As opposed to only 0.131 > n^h"^ being imposed.) 

If CDMS does not make a positive detection by the end of 2007, there will be implications 
for heavy MSSM Higgs searches at the Tevatron and LHC. Similarly, the lack of a detection 
of A/H + X — s> r+r^ + X at the Tevatron will impact the prospects for future direct dark 
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FIG. 6: The same as in Fig. [5l but showing only those supersymmetric models which are not within 
the 2007 projected reach of CDMS (left frames) or which are not within the projected reach of the 
Tevatron (after 4 fb~^) for the process A/H + X — > t~^t~ +X (right frames). As in Fig. [5l the upper 
frames show those models which do not overproduce the thermal abundance of neutralino dark matter 
(0.131 > n^/i^) while the lower frames also require that the models shown do not predict a smaller 
thermal abundance of neutralino dark matter than measured by WMAP (0.131 > ^x^'^ ^ 0.089) |:2^.) 

matter searches. In Fig. [6l we plot (in a fashion similar to Fig. [5]) the supersymmetric models 
found in our scan which are not within the 2007 projected reach of CDMS (left frames) or are 
not within the reach of the Tevatron (after 4 fb~^) for the process A/H + X —>■ t^t" + X 
(right frames). If CDMS does not make a positive detection by the end of 2007, only a small 
fraction of the remaining models are within the reach of the Tevatron through the channel 
A/H + X r^r" + X. These models are those in which the lightest neutralino has a very 
small higgsino fraction (large |/i|). On the other hand, if the Tevatron does not observe the 
process A/H + X t^t~ +X, the reach of direct detection experiments can still be promising. 
Whether heavy MSSM Higgs bosons are observed at the LHC will also have implications for 



15 



10" 
10" 

10" 

10" 

10" 

10" 

10"' 
10-' 







1 

jFmi - 

■ -rt5ilS 




^ " ^ cow's 20°^- 
























100 


1000 




1000 



10" 
10" 
10" 
10" 
10" 
10" 

10-' 
10-' 



1 1 1 1 1 1 1 


1 


"5«*^« ^^''^ ^"if^' 


' . _>1^— -"^ ^ 








_ _ , - ■ f v^-v- ; 
















^x '»A*^" 




, , , , 1 





10- 
10- 




. . . I- 1 • ■■. ■ L. , . . . 



100 



1000 



100 



1000 



FIG. 7: The same as in the right frames of Figs. [5] and [U but showing only those super symmetric 
models which are (left frames) and are not (right frames) within the projected reach of the LHC 
(after 30 fb^^). In each frame, those models within the reach of the Tevatron are not shown. As in 
Figs. El and m the upper frames show those models which do not overproduce the thermal abundance 
of neutralino dark matter (0.131 > O^^/i^) while the lower frames also require that the models shown 
do not predict a smaller thermal abundance of neutralino dark matter than measured by WMAP 
(0.131 > n^h'^ > 0.089) H.) 



the prospects of future direct detection experiments. In Fig. [7] we plot, in the ~ "^x plane, 
those models found by our scan which are (left frames) and are not (right frames) within the 
reach of the LHC heavy neutral MSSM Higgs search (with 30 fb~^). Once again, the upper 
figures require 0.131 > Q^h'^ while the lower frames require 0.131 > ^^^"^ ^ 0.089. From the 
left frames, we conclude that if a positive detection is made at the LHC, then the prospects 
for directly detecting neutralinos with CDMS or the early phases of Super-CDMS will be very 
promising. 

The combination of these figures demonstrates that searches at the LHC and direct dark 
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matter searches can be highly complementary. A wide range of supersymmetric models exist 
in which one, both or neither of the LHC and direct detection experiments will be successful 
in observing heavy MSSM Higgs bosons, or neutralino dark matter, respectively. 

Note that here we have only considered 30 fb~^ of data for the LHC. Further improvements 
on the reach of the LHC with greater luminosity are to be expected. 



VI. UNCERTAINTIES, CAVEATS AND LIMITATIONS 

Throughout our study, we have adopted a number of assumptions. In this section, we 
summarize these, and discuss how our conclusions might be affected by relaxing them. 



A. Astrophysical Uncertainties 

Throughout our study, we have adopted a standard local dark matter density (p = 0.3 
GeV/cm^) and velocity distribution (an isothermal Maxwell-Boltman distribution). If the 
actual local density or velocity distribution of dark matter is substantially different, then clearly 
our conclusions may be modified. 

The local density and velocity distribution of dark matter can be inferred by studying the 
rotation curves of our galaxy. Different groups have come to somewhat different conclusions 
regarding these observations: Bahcall et al. find a best-fit value of p = 0.34GeV/cm^ j29| . 
Caldwell and Ostriker find p = 0.23GeV/cm^ jsO] while Gates, Gyuk and Turner find p = 



0.34 — 0.73GeV/cm [3l|. Bergstrom, UUio and Buckley [32| find that the observations are 
consistent with local dark matter densities in the range of about 0.2 — 0.8GeV/cm^. The 
uncertainties in the local dark matter velocity distribution are less important in estimating the 
rates in direct dark matter detection experiments [33|. 

Rotation curves, however, only constrain the dark matter density as averaged over scales 
larger than a kiloparsec or so. In contrast, the solar system moves a distance of ~10~'^ parsecs 
relative to the dark matter halo each year. If dark matter is distributed in an inhomogeneous 
way over milliparsec scales {ie. as a collection of dense clumps and voids), then the density 
along the path of the Earth, as seen by direct detection experiments, could be much larger or 
smaller than is inferred by the rotational dynamics of our galaxy. This is not anticipated to 
be a problem, however. The vast majority of dark matter in the inner regions of our galaxy 
have been in place for ~10^° years; ample time for the destruction of clumps through tidal 
interactions. Using high-resolution simulations, Helmi, White and Springel find that the dark 
matter in the solar neighborhood is likely to consist of a superposition of hundreds of thousands 
of dark matter streams, collectively representing a very smooth and homogeneous distribution 



34 



If the uncertainties in the local halo density, the velocity distribution, and the parameters 
/jf '"^ (from Eq. El are taken into account, the reach of direct dark matter searches are not a 
single line, but instead should be thought of as a band whose width can be a factor of 3 to 5, 
depending on which estimates are adopted for the range of acceptable local densities. 

Finally, regarding the density of dark matter, thoughout this study we have assumed that 
all of our universe's dark matter consists of neutralinos. If neutralinos make up only a fraction 
of the local dark matter density, then the rates in direct detection experiments will be reduced 
accordingly. 
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B. Assumptions in the Supersymmetric Model 



As we do not know through what mechanism or mechanisms supersymmetry is broken, 
there is a extremely vast range of characterstics the supersymmetric spectrum might possess. 
To attempt to study the entire range of parameter space within the MSSM is not generally 
a tractable approach. If extensions of supersymmetry beyond the MSSM are considered, this 
becomes only more difficult. For this reason, we have adopted a number of assumptions regard- 
ing the nature of supersymmetry. In particular, we have limited ourselves to the MSSM with 
no CP-violating phases, and have adopted the GUT-relationship between the gaugino masses 
{Ml/ gi = M2/g2 = M^/g^). We have also adopted five common soft SUSY breaking sfermion 
mass parameters {mQ., mu., rriD., m^. and m^;.) at the low energy scale. 

Adoping the GUT relationship between the gaugino masses guarantees that the wino-fraction 
of the lightest neutralino will be quite small (a few percent or less). From Eqns. [TOl and [TT| 
we see that the wino fraction of the lightest neutralino, if large, can play an important role in 
determining its elastic scattering cross section with nucleons. In the case that M2 -C Mi, in 
addition to the lightest neutralino being wino-like, the lightest chargino is only slightly heavier 
than the LSP. In such a scenario, neutralino-neutralino annihilations and neutralino-chargino 
coannihilations are too efficient to produce the measured abundance of thermal neutralinos 
(unless ^ M2 > 3 TeV). If all or most of the dark matter consists of neutralinos in 
such a model, non-thermal mechanisms must be relied upon. Non-thermally generated wino 
dark matter can appear naturally in models of Anomaly Mediated Supersymmetry Breaking 
(AMSB) for example. 

We have also not considered the possibility of significant CP-violating phases in the MSSM. 
The presence of such phases can impact the neutralino-nucleon elastic scattering cross section, 
the prospects for heavy Higgs discovery at colliders, and the neutralino relic abundance through 



their effect on the sparticle and Higgs boson masses and couplings [36 . 

Lastly, if we consider supersymmetric models beyond the MSSM, a very broad range of pos- 
sibilities become available. Neutralino dark matter has been studied in the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM) and the near-Minimal Supersymmetric Standard 
Model (nMSSM) by a number of authors [37|. The neutralino spectrum in these extensions 
contain, in addition to the four neutralinos found in the MSSM, a singlino (the superpartner 
of an additional Higgs singlet), leading to modified couplings and masses for the lightest neu- 
tralino. The CP-even singlet Higgs in these extensions could also mediate neutralino-nucleon 
elastic scattering, for example. In other extensions of the MSSM, many diverse neutralino 
phenomenologies can be found {38i] . 



VII. DISCUSSION AND CONCLUSIONS 

In this article, we have studied the interplay between collider searches for MSSM Higgs 
bosons and direct dark matter experiments. In particular, the prospects for heavy MSSM Higgs 
searches at the Tevatron and LHC and for direct dark matter searches are each most promising 
in the case of large values of tan/3 and small values of m^. There is, therefore, an interesting 
relationship between these two classes of experiments. Our findings can be summarized by the 
following: 

• If neutralinos are detected by direct dark matter experiments in the near future, then the 
prospects for MSSM heavy Higgs searches at the Tevatron and LHC will be significantly 
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enhanced. In this case, however, it remains possible that the large neutralino-nucleon 
elastic scattering cross section could result largely from the composition of the lightest 
neutralino (a mixed bino-higgsino, for example) rather than from itla and tan (3 being 
light and large, respectively. 

• If the Tevatron or LHC observes heavy MSSM Higgs bosons, then the prospects for the 
direct detection of neutralinos will be excellent. Models which are observable at the 
Tevatron in these channels (with 4 fb~^ luminosity) typically predict neutralino-nucleon 
cross sections in the range of 10~^ to 10~^ pb. Models which are observable at the LHC in 
these channels (with 30 fb~^ luminosity) have a wide range of possible elastic scattering 
cross sections, roughly 10"^ to 10^^^ pb, depending on the composition of the lightest 
neutralino. 

• If no WIMPs are detected by direct dark matter experiments in the near future (by CDMS 
by the end of 2007, for example), then MSSM heavy Higgs searches at the Tevatron 
are expected to be potentially successful only if the lightest neutralino is very bino-like 
(|/i| > 800 GeV).^ The prospects for heavy MSSM Higgs discovery at the LHC are 
promising even if no signal is observed by direct detection experiments. 

• If the Tevatron and/or LHC do not observe heavy MSSM Higgs bosons, then the prospects 
for the direct detection of neutralinos will be reduced, although many models will remain 
viable which are within the reach of near future and planned direct dark matter searches. 

In this work we have explored one aspect of the interplay between collider experiments and 
astrophysical observations, which relies on testing different sectors of the MSSM relevant for the 
understanding of dark matter. We would like to emphasize the importance of using both collider 
and direct detection experiments to study supersymmetry. If the neutralino is discovered at 
the Tevatron or the LHC, it will be through missing energy signatures in combination with jets 
and/or leptons. The LSP candidate lifetime would be constrained to microsecond time scales 
or longer, but the possibility that this particle is only metastable, and thus not abundant in the 
universe today, would remain open. This could occur as a result of R-parity violating couplings, 
for example, or through neutralino decays to a lighter gravitino or axino [s^. Until the same 
particles are observed in terrestial direct detection experiments, it will be impossible to know 
that they are stable over cosmological time scales. 

The LHC experiments will search for evidence of dark matter particles in events with large 
missing energy plus multiple jets and/or leptons, as expected in the cascades from heavy colored 
particles like gluinos and squarks, for example. If the new colored particles are within the reach 
of the LHC, namely if their masses are below a few TeV, then it is likely that the LHC will find 
evidence for dark matter particles. In the simplest models analyzed, this corresponds to WIMP 
masses up to a few hundred GeV. In scenarios in which these colored particles are too heavy 
to be produced at the LHC, the direct production of other new, weakly interacting particles, 
which ultimately decay into the dark matter candidate, remains as a possible search channel. 
In these cases, multi-lepton signals plus missing energy are the most robust option against the 



^ Although Figs. [3] and 2] show a value of > 1200 GeV rather than 800 GeV, this is the result only for the 
special case of A/i ~ 100 GeV, as is used in those figures. The value derived from a general scan in the Mi-/x 
plane is l^i] > 800 GeV [lO|. 
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copious QCD backgrounds. In many case studies it appears that the properties and couphngs 
of the new particles can be measured with, at best, only modest precision, and hence it will be 
difficult to establish the identity of the dark matter particle from LHC data alone. Direct dark 
matter detection experiments would open another window into the nature and composition of 
neutralinos, and hence play an essential and complementary role in the quantitative study of 
super symmetric dark matter (4o| . 

In the more distant future, experiments at the prospective International Linear Collider 
(ILC) would measure many of the properties of supersymmetric particles much more accu- 



rately than can be done at the LHC [41|. In the case that the non-strongly interacting particles 
associated with the dark matter are within the kinematic reach of the ILC, these particles 
can be produced and their masses and couplings measured with high precision. The avail- 
ability of polarized beams and the capability to make precise measurements of cross sections 
is particularly useful. The ILC can also give direct information on particle masses from the 
kinematic distributions of decay products and from the measurement of excitation curves at 
center-of-mass energies near threshold. In the case of heavy new particles that are weakly 
interacting, the above measurements can constrain the heavy masses and couplings relevant to 
the computation of the dark matter relic density. Ultimately, information from the ILC can 
determine the relevant supersymmetric parameters with sufficient accuracy such that the relic 
density can be computed to the percent level and compared with cosmological measurements 
of similar precision. Given an understanding of the properties of the neutralino at this level, 
measurements of rates at direct (and indirect) dark matter detection experiments will allow 



one to infer the local distribution of dark matter |42l . l43l |. and even constrain the expansion 



history of our universe since the time of neutralino freeze-out. 



Acknowledgments 

DH and AV are supported by the Department of Energy and by NASA grant NAG5-10842. 
MC is supported by the US Department of Energy grant DE-AC02-76CHO3000. 



[1] D. S. Akerib et al. [CDMS Collaboration], Phys. Rev. Lett. 96, 011302 (2006) 
|arXiv:astro-pli/0509 259|; D. S. Akerib et al. [CDMS Collaboration], Phys. Rev. D 73, 011102 
(2006) ; arXiv:astro-ph /0509269 1 . 

[2] G. J. Alner et al. [UK Dark Matter Collaboration], Astropart. Phys. 23, 444 (2005). 

[3] V. Sanglard et al. [The EDELWEISS Collaboration], Phys. Rev. D 71, 122002 (2005) 
|arXiv:astro-ph/0503265|. 

[4] G. Angloher et al., Astropart. Phys. 23, 325 (2005) [ arXiv:astro-ph/ 0408006| . 

[5] R. Brunetti et al, New Astron. Rev. 49, 265 (2005) arXiv:astro-ph/0405342] ; Preliminary lim- 
its using liquid Argon were shown at the April 2006, P5 meeting at Fermilab by Elena Aprile 
(http://www.fnal.gov/directorate/program_planning/P5/P5_Apr2006/Talks/Aprile.pdf). 

[6] G. Bertone, D. Hooper and J. Silk, Phys. Rept. 405, 279 (2005) | arXiv:hep-ph /0404175 1 . 

[7] D. Bortoletto [CDF and DO Collaborations], PoS HEP2005, 347 (2006); A. Canepa [CDF Col- 
laboration], arXiv:hep-ex/0603032. 

[8] A. Anastassov [CDF and DO Collaborations], PoS HEP2005, 326 (2006); 



20 



[9] V. Abazov [DO Collaboration], |arXiv:hep-ex /0604029 , 
[10] M. Carena, D. Hooper and P. Skands, Phys. Rev. Lett. 97, 051801 (2006) [arXiv:liep-ph/0 603180] . 
[11] M. Carena, S. Heinemeyer, C. E. M. Wagner and G. Weiglein, arXiv:hep-ph/9912223. 
[12] S. Heinemeyer, W. Hollik and G. Weiglein, Comput. Phys. Commun. 124, 76 (2000) 
I arXiv:hep-p li/98123 20] . 

[13] J. S. Lee, A. Pilaftsis, M. Carena, S. Y. Choi, M. Drees, J. R. Ellis and C. E. M. Wagner, Comput. 

Phys. Commun. 156, 283 (2004) |arXiv:hep-piV0307377] . 
[14] S. P. Martin, arXiv:hep-ph/9709356', 

[15] M. Carena, R. L. Culbertson, S. Eno, H. J. Frisch and S. Mrenna, Rev. Mod. Phys. 71, 937 

(1999) |arXiv:hep-ex/9712022l . 
[16] M. Carena and H. E. Haber, Prog. Part. Nucl. Phys. 50, 63 (2003) [arXiv:hep-ph/02 08209]. 
[17] A. Bottino, F. Donato, N. Fornengo and S. Scopel, Astropart. Phys. 18, 205 (2002) 

[arXiv:hep-ph/0 1 11229] ; Astropart. Phys. 13, 215 (2000) [arXiv:hep-ph/9909228]; J. R. El- 

hs, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Rev. D 71, 095007 (2005) 

|arXiv:hep-ph/0502001|. 

[18] G. B. Gelmini, P. Gondolo and E. Roulet, Nucl. Phys. B 351, 623 (1991); M. Srednicki and 
R. Watkins, Phys. Lett. B 225, 140 (1989); M. Drees and M. Nojiri, Phys. Rev. D 48, 3483 
(1993) |arXiv:hep-ph/9307208|; M. Drees and M. M. Nojiri, Phys. Rev. D 47, 4226 (1993) 
arXiv:hep-ph/9210272]; J. R. Ellis, A. Ferstl and K. A. Olive, Phys. Lett. B 481, (2000) 304, 



arXiv:hep-ph/0001005 



[19] A. Abulencia et al. [CDF Collaboration], Phys. Rev. Lett. 96, 011802 (2006) 
|arXiv:hep-ex/0508051 1; A. Abulencia et al. [CDF Collaboration], Phys. Rev. Lett. 96, 
042003 (2006) [arXiv:hep-ex/0510065 1. 

[20] V. M. Abazov et al. [DO Collaboration], |arXiv:hep-ex/0605009; V. M. Abazov et al. 



[DO Collaboration], Phys. Rev. Lett. 95, 151801 (2005) [arXiv:hep-ex/0504018] ; www- 
dO.fnal.gov/Run2Physics/WWW/results/prehm/HIGGS/H24/H24.pdf. 
[21] M. Carena, S. Heinemeyer, C. E. M. Wagner and G. Weiglein, Eur. Phys. J. C 45, 797 (2006) 
[arXiv:hep-ph/0511023|. 

[22] G. W. Bennett et al. [Muon g-2 Collaboration], Phys. Rev. Lett. 92, 161802 (2004) 
I arXiv: hep-ex/0401008 1 ; R. R. Akhmetshin et al. [CMD-2 Collaboration], Phys. Lett. B 578, 
285 (2004) pirX"iv:hep-ex/0308008l ; M. Davier, S. Eidelman, A. Hocker and Z. Zhang, Eur. Phys. 
J. C 31, 503 (2003) [arXiv:hep-ph/0308213|; F. Jegerlehner, Nucl. Phys. Proc. Suppl. 126, 325 



(2004) [arXiv:hep-ph/0310234 


; K. Hagiwara, A. 


Rev. D 69, 093003 (2004) 


arXiv:hep-ph/0312250 



[23] T. Moroi, Phys. Rev. D 53, 6565 (1996) [Erratum-ibid. D 56, 4424 (1997)] 
|arXiv:hep-ph/9512396| ; M. Carena, G. F. Giudice and C. E. M. Wagner, Phys. Lett. B 
390, 234 (1997) |arXiv:hep-ph/9610233|; L. L. Everett, G. L. Kane, S. Rigolin and L. T. Wang, 
Phys. Rev. Lett. 86, 3484 (2001) [arXiv :hep-ph /01 02145| ; E. A. Baltz and P. Gondolo, Phys. 
Rev. D 67, 063503 (2003) [arXiv:astro-ph/0207673|. 

[24] M. S. Alam et al. [CLEO Collaboration], Phys. Rev. Lett. 74, 2885 (1995); K. Abe et al. [Belle 
Collaboration], |arXi v:hep-ex/0l 07065 ; L. Lista [BABAR Collaboration] , |arXiv:he p-ex/01100To| 
B. Aubert et al. [BABAR Collaboration], arXiv:hep-ex/0207074; 

[25] R. Barbieri and G. F. Giudice, Phys. Lett. B 309, 86 (1993) [arXiv:hep-ph/9303270]; S. Bertohni, 
F. Borzumati, A. Masiero and G. Ridolfi, Nucl. Phys. B 353, 591 (1991); M. Carena, D. Garcia, 
U. Nierste and C. E. M. Wagner, Phys. Lett. B 499, 141 (2001) [arXiv:hep-ph/0010003l; G. De- 



21 



grassi, P. Gambino and G. F. Giudice, JHEP 0012, 009 (2000) |arXiv:hep-ph /0009337|; H. Baer, 

M. Brhlik, D. Castano and X. Tata, Phys. Rev. D 58, 015007 (1998) |a?Xiv:hep-ph/9712305]; 

M. Carena, A. Menon, R. Noriega-Papaqui, A. Szynkman and C. E. M. Wagner, Phys. Rev. D 

74, 015009 (2006) |arXiv:he p-ph/0603106| . 
[26] CDF exotics public webpage: Ihttp: / / www-cdf.fnal.gov/physics/exotic/r2a/20050519.mssm_ htt / projections. h 
[27] P. Gondolo, J. Edsjo, P. Ulho, L. Bergstrom, M. Schelke and E. A. Baltz, New Astron. Rev. 49, 

149 (2005); P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schelke and E. A. Baltz, JCAP 

0407, 008 (2004) |arXiv:astro-ph/0406204]. 
[28] D. N. Spergel et al. [WMAP Cohaboration] , arXiv:astro-ph/0603 449l 
[29] J. N. Bahcah, M. Schmidt and R. M. Soneira, Astrophys. J. 265 (1983) 730. 
[30] R. R. Caldwell and J. P. Ostriker, Astrophys. J. 251 (1981) 61. 



[31] E. I. Gates, G. Gyuk and M. S. Turner, Astrophys. J. 449, L123 (1995) [ast ro-ph/95 05039]. 
[32] L. Bergstrom, P. Ulho and J. H. Buckley, Astropart. Phys. 9 (1998) 137 [arXiv:astro-ph/9 712318j . 
[33] M. Kamionkowski and A. Kinkhabwala, Phys. Rev. D 57, 3256 (1998) [ hep-ph/9710337 ]; 
[34] A. Helmi, S. D. M. White and V. Springel, Phys. Rev. D 66, 063502 (2002) 
iarXiv:astro-ph/0201289 | . 

[35] T. Moroi and L. Randall, Nucl. Phys. B 570, 455 (2000) [arXiv:hep-ph/9906527]; P. Ullio, JHEP 
0106, 053 (2001) [arXiv:hep-ph/ 0105052 |; D. Hooper and L. T. Wang, Phys. Rev. D 69, 035001 
(2004) [a rXiv:hep-ph/0309036] . 

[36] C. Balazs, M. Carena and C. E. M. Wagner, Phys. Rev. D 70, 015007 (2004) 
|arXiv:hep-ph/0403224] ; C. Balazs, M. Carena, A. Menon, D. E. Morrissey and C. E. M. Wagner, 
Phys. Rev. D 71, 075002 (2005) | arXiv:hep-ph /0412264 1 ; G. Belanger, F. Boudjema, S. Kraml, 
A. Pukhov and A. Semenov, Phys. Rev. D 73, 115007 (2006) [arXiv:hep-ph/0604150]; T. Ni- 
hei and M. Sasagawa, Phys. Rev. D 70, 055011 (2004) [Erratum-ibid. D 70, 079901 (2004)] 
|arXiv:hep-ph/0404100 1 . 

[37] D. G. Cerdeno, C. Hugonie, D. E. Lopez-Foghani, C. Munoz and A. M. Teixeira, JHEP 0412, 
048 (2004) |arXiv:hep-ph/0408102|; G. Belanger, F. Boudjema, C. Hugonie, A. Pukhov and 

A. Semenov, JCAP 0509, 001 (2005) [arXiv:hep-p h/05051 42|; J. F. Gunion, D. Hooper and 

B. McEhath, Phys. Rev. D 73, 015011 (2006) [arXi v:hep-ph/0509024] ; V. Barger, P. Langacker 
and G. Shaughnessy, arXiv:hep-ph/0609068; A. Menon, D. E. Morrissey and C. E. M. Wagner, 
Phys. Rev. D 70, 035005 (2004) |arXiv:hep-ph/0404184 |. 



[38] V. Barger, P. Langacker and H. S. Lee, Phys. Lett. B 630, 85 (2005) [arXiv:hep-ph/0508027]; 
V. Barger, C. Kao, P. Langacker and H. S. Lee, Phys. Lett. B 600, 104 (2004) 
[ arXiv:hep-ph/04081 20| . 

[39] J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. Lett. 91, 011302 (2003) 
|arXiv:hep-ph/0302215] ; J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. D 68, 063504 

(2003) iarXi v:hep-ph /0306024 1 ; J. L. Feng, S. Su and F. Takayama, Phys. Rev. D 70, 075019 

(2004) [arXiv:hep-ph/0404231j. 



[40] D. Hooper and A. M. Taylor, |arXiv:hep-ph/0607086 



[41] G. Weiglein et al. [LHC/LC Study Group], Phys. Rept. 426, 47 (2006) |arXiv:hep-ph/0410364] . 
[42] E. A. Baltz, M. Battagha, M. E. Peskin and T. Wizansky, |arXiv:hep-ph / 0602187t 
[43] M. Carena and A. Freitas, arXiv:hep-ph/0608255i M. Carena, A. Finch, A. Freitas, C. Milstene, 
H. Nowak and A. Sopczak, Phys. Rev. D 72, 115008 (2005) [arXiv:hep-ph/0508152j . 



22 



